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The linear reversal mechanism in FePt grains ranging from 2.316 nm to 5.404 nm has been simulated using 
atomistic spin dynamics, parametrized from ab-initio calculations. The Curie temperature and the critical 
temperature (T*), at which the linear reversal mechanism occurs, are observed to decrease with system size 
whilst the temperature window T* < T < Tq increases. The reversal paths close to the Curie temperature 
have been calculated, showing that for decreasing system size the reversal path becomes more elliptic at lower 
temperatures, consistent with the decrease in the Curie temperature arising from finite size effects. Calcula¬ 
tions of the minimum pulse duration show faster switching in small grains and is qualitatively described by 
the Landau-Lifshitz-Bloch equation with finite size atomistic parameterization, which suggests that multiscale 
modeling of FePt down to a grain size of « 3.5 nm is possible. 


The properties of Llo FePt are of paramount impor¬ 
tance for the next generation of high density magnetic 
recording technology, in particular Heat Assisted Mag¬ 
netic Recording (HAMR). Llo FePt exhibits a very large 
uniaxial anisotropy of K u w 7 x 10 6 Jm -3 which allows 
possible stable grain sizes down to a limit of around 3 
nm. 1 To record data the medium is heated close to or 
above the Curie temperature where the anisotropy is sig¬ 
nificantly reduced, therefore an understanding of the dy¬ 
namical behavior of the nanometer FePt grains at ele¬ 
vated temperatures is essential. In particular, at tem¬ 
peratures close to Tc it has been shown that the mag¬ 
netization reverses by a linear longitudinal route rather 
than a precessional one. 2,3 Using the Landau-Lifshitz- 
Bloch model Kazantseva et al. 2 predicted that the linear 
reversal regime occurs at a critical temperature T* and 
derived a expression for the reversal time in this regime. 
Following this Barker et al. 3 observed the linear rever¬ 
sal mechanism using an atomistic model of FePt which 
agreed well with the LLB modeling. However the ques¬ 
tion of linear reversal in finite nanometer grains and the 
size dependence of the switching properties remains open. 

Here an atomistic model is used to investigate the 
properties of nanometer grains from 2.316 nm to 5.404 
nm at temperatures close to the Curie temperature. Us¬ 
ing Langevin dynamics methods we investigate the static 
and dynamic properties as a function of the grain size. 
We calculate the equilibrium magnetization and suscepti¬ 
bilities for each system size which allow us to parametrize 
the LLB model. Using this parametrization we can then 
compare the reversal path ellipticity and times calculated 
using the atomistic model to those predicted by the LLB. 

The properties of nanometer FePt particles have been 
investigated using ab-initio methods by Chepulskii and 
Butler. 4 Their calculations show that the magneto¬ 
crystalline anisotropy and magnetic moment varies sig¬ 
nificantly with system size and shape. However, the sizes 
considered are smaller than 1 nm so here we use the pa¬ 
rameterization of bulk materials given by Mryasov et al. 5 
The finite size effects in FePt have been investigated us¬ 
ing a similar model by Lyberatos et al. 6 and Hovorka 


et al? Lyberatos et al. have further investigated the 
grain size effects on the HAMR write process showing 
that smaller grain sizes require larger write fields. 8 

Atomistic spin models are increasingly used for dy¬ 
namic magnetization calculations (Evans et al. 9 ). The 
ab-intio simulations of Mryasov et al. show that the Llo 
FePt spin Hamiltonian can be written in the following 
form, involving only Fe degrees of freedom with parame¬ 
ters mediated by the Pt: 
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The first term is the mediated exchange interaction, the 
second term a uniaxial anisotropy and the third a two- 
ion anisotropy which represents an anisotropic exchange 
interaction. The sum j is over the atoms within the 
interaction range, which extends for approximately 5 
unit cells. It is important to note that our model in¬ 
cludes a significant finite size effect in that, as noted 
by Nowak et al. 10 , there is a loss of anisotropy at the 
surface due to the reduction of coordination which re¬ 
sults in a reduction of the 2-ion anisotropy contribution. 
These finite size effects mean the total anisotropy varies 
from K u = 8.3 x 10 6 Jm~ 3 for a 4.632 nm particle to 
K u = 7.3 x 10 6 Jm -3 for a 2.316 nm particle. The dipo¬ 
lar held was not included in the calculations as at tem¬ 
peratures close to Tc it has negligible effect and is not 
important in the linear reversal mechanism. 

The time evolution of the magnetic properties is cal¬ 
culated using the Langevin dynamic approach. The LLG 
equation for an atomic spin is: 

f = arku s ‘ x(H - +ls - xH - > (2) 

Where 7 = 1.76 x 10 n s _1 T _1 is the gyromagnetic ra¬ 
tio, n s = 3.23 )jlb is the saturation magnetic moment of 
the combined Fe and Pt atoms, A = 0.1 is the atomistic 
damping/coupling parameter and H; is the effective held 
acting on the spin. For these parameters the saturation 
magnetization is M s = 1.075 x 10 6 JT _1 m -3 . Temper¬ 
ature is introduced via a thermal noise term which is 
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FIG. 1. The temperature dependence of the magnetization 
for particle sizes of 2.316 nm to 5.404 nm. The dashed black 
line shows the magnetization for a bulk system from ref. 14 
for comparison. As the system size decreases the magnetiza¬ 
tion below the Curie temperature decreases while above it the 
magnetization increases. The inset shows the Curie tempera¬ 
tures which have been extracted from fitting and T* extracted 
from figure 2. The solid lines are equation 5 fitted to find the 
critical exponent. 


added to the effective field, = —dW/dSi + The 
thermal noise is a Gaussian stochastic process with the 
following mean and variance. 11 
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where i,j refer to atom numbers, a, (3 Cartesian coordi¬ 
nates and t, t' the time. The resulting stochastic LLG 
equation is numerically solved using the Semi-Implicit 
integration scheme 12,13 with a time step At = 10 -16 s 
which has been tested to show stability. 

We first investigate the temperature dependence of the 
magnetization and the apparent Curie temperature. The 
results are shown in figure 1 for a selection of system sizes. 
Clearly the smaller systems have a reduced equilibrium 
magnetization at a given temperature, as expected as 
a result of reduced co-ordination at the surfaces. The 
apparent Curie temperature can be extracted from fitting 
m{T) = (1 — T/Tc{L)Y to the magnetization curves. 6 
The extracted Tc values are reduced significantly from 
the bulk and described well by finite system size scaling 
theory: 7 
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The inset in figure 1 shows the fitting to obtain TJ2? = 665 
K, v = 0.85699, do = 0.4039 nm which agree well with 
refs. 7 and 6 for the long range FePt Hamiltonian which 
we utilize. 


FIG. 2. The ratio of the parallel and perpendicular suscepti¬ 
bilities reach a peak at different temperatures consistent with 
the Curie temperature extracted from the magnetization. The 
solid lines show functions that have been fitted the susceptibil¬ 
ities and the horizontal solid line gives X||/X-L = 1/2 defining 
the transition to linear reversal. 


The switching time at elevated temperatures is 
strongly affected by the reversal mechanism, which is de¬ 
termined by the ratio of longitudinal and transverse sus¬ 
ceptibilities, X|| , X-L- Consequently we next investigate 
the ratio of susceptibilities and make a comparison with 
bulk properties using an expression for the free energy, 
in the absence of an externally applied field 15,16 : 
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where m is the reduced magnetization, m e the equilib¬ 
rium magnetization and X||,X_L are respectively the re¬ 
duced parallel and perpendicular susceptibility. We note 
that eqn. 6 gives the free energy governing the behavior 
of the LLB equation derived in ref. 15, which allows us to 
connect the current calculations to the LLB equation as 
will be shown later. From the free energy we find the fol¬ 
lowing expression for the orientation dependence of the 
magnetization length, m(9) for a bulk material: 
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Eqn. 7 gives an important relationship between the ellip- 
ticity parameter and the characteristic parameters of the 
LLB equation. Specifically, the ratio of the longitudinal 
and perpendicular susceptibility determines the clliptic- 
ity of the reversal path. Eqn. 7 predicts a critical tem¬ 
perature, T*, corresponding to X\\(T*)/xj_(T*) = 1/2, 
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at which point the perpendicular magnetization vanishes. 
T* is the critical temperature above which the magnetiza¬ 
tion will reverse through the linear switching mechanism. 
By calculating the equilibrium magnetization and suscep¬ 
tibilities for finite nanometer systems equations 6 and 7 
can be parametrized. The longitudinal and transverse 
susceptibilities have been calculated as Xz and \x + Xy 
respectively. The ratio x\\(T)/x±(T) is shown in fig 2. 
The points show the susceptibilities from the atomistic 
simulations while the solid lines show a polynomial func¬ 
tion that has been fitted to inverse susceptibilities in¬ 
dividually in a similar manner to ref. 14. The signifi¬ 
cant noise in the small grains arises from the thermal 
switching of the magnetization during the averaging of 
rn z leading to ergodicity breaking. 1 ' As the grain size 
decreases the critical temperature T* decreases and the 
apparent temperature window T* < T < Tq of the linear 
reversal region widens. T* is shown in the inset of fig 1 
alongside Tc which highlights the widening temperature 
window. Although the free energy expression given in 
eqn. 6 is strictly valid only for bulk systems, the tem¬ 
perature variation of X|| (T)/x_l(T) for all system sizes is 
consistent with the behavior shown by the reversal paths 
and the ellipticity. 

Figure 3.(a)-(d) shows the reversal paths for system 
size 2.314 nm to 4.632 mn using a 1 T field to investi¬ 
gate the temperature dependence of the reversal path. In 
the 4.632 nm particle the magnetization does not reverse 
within the simulation time limit of 400 ps for 540 K and 
560 K but then there is a sharp transition to linear rever¬ 
sal, characterized by vanishing transverse magnetization, 
at 640K. The smaller systems reverse at lower tempera¬ 
tures but there is a less critical transition to the linear 
reversal path. 

To further investigate the effect of the system size on 
the transition to the linear regime, we calculate the el¬ 
lipticity of the reversal paths, defined as 
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The ellipticity is shown in figure 3.(e) for the reversal 
paths shown in fig. 3.(a)-(d). As discussed previously 
the transition to the linear regime (where e « 1) oc¬ 
curs within a smaller temperature window for the larger 
systems. Also since the Curie temperature decreases for 
smaller system sizes the window is centered at lower tem¬ 
peratures. The lines shown in fig 3.(e) show the elliptic¬ 
ity from the LLB model which is cllb = ^X\\ (T)/x±(T), 
this shows that the LLB model reasonably predicts the 
ellipticity but significant deviations occur for the smaller 
particles. 

In terms of HAMR, the temperature dependence of the 
relaxation time is an important factor. It has been shown 
by Evans et al. 18 that the achievable recording density 
in HAMR is set by the equilibrium magnetization at the 
temperature at which the magnetization freezes. Con¬ 
sequently, it is vital to the HAMR process for the mag¬ 
netization to be as close as possible to the equilibrium 


value at a given temperature, emphasizing the impor¬ 
tance of the linear reversal mechanism and its associated 
fast switching time governed by longitudinal rather than 
transverse relaxation. As a measure of the relaxation 
time we calculate the minimal pulse duration 3 ; this is 
the time taken for rn z to first pass the m z = 0 plane 
starting from an initial fully magnetized state. Figure 4 
shows the minimal pulse duration calculated with a 1 T 
reversing field relevant to the HAMR process. The results 
show that there is a significant reduction in the relaxation 
time in small grains which arises from a convolution of 
two factors, the decrease of intrinsic properties and the 
transition from circular to linear reversal. The intrinsic 
parameters K and M s both decrease with temperature 
and depend strongly on the system size. Since a large 
part of the anisotropy arises from the 2-ion interaction 
there is a strong decrease with grain size and thus the re¬ 
laxation time predictable from the Arrhenius-Neel law. 19 
In the elliptical and linear regimes the free energy barrier 
is reduced significantly relative to the coherent reversal 
mechanism, leading to a further reduction of switching 
time with temperature over and above that predicted by 
the Arrhenius-Neel law. We note that both factors are es¬ 
sentially thermodynamic and emerge naturally from the 
atomistic model. 

We compare these results to the expression for the min¬ 
imum pulse duration derived from the LLB equation by 
Kazantseva et al. 2 The minimum pulse duration depends 
on the equilibrium magnetization and susceptibilities so 
using the values already presented it can be computed for 
each specific grain size, which are shown as solid lines in 
figure 4. There is a good agreement but since the free en¬ 
ergy assumes an infinite system the problems arise at size 
dependent Tc so only the section above is shown. The 
analytic expression with finite input parameters shows 
qualitative agreement with the atomistic reversal times 
showing that for moderate system sizes the LLB equa¬ 
tion is applicable using parameters determined for the 
specific system size. The reversal times for the 6.948 nm 
particle are similar to the bulk LLB curve. For the small¬ 
est system sizes the validity is questionable since it is far 
from the bulk regime in which the LLB is derived. How¬ 
ever, for the sizes in excess of 3.5 nm potentially usable 
as HAMR media, the parameterization seems a practical 
proposition. 

In summary we have investigated the effects of finite 
size on the linear reversal regime of FePt which is of 
technological importance for HAMR. Atomistic calcu¬ 
lations of the temperature dependence of the magneti¬ 
zation shows that Tc decreases with decreasing system 
size as does the critical temperature for the linear re¬ 
versal regime, T* determined from the susceptibilities. 
The transition to linear reversal also takes place at a 
lower temperature, with the temperature range of lin¬ 
ear reversal increasing with decreasing size. The reversal 
characteristics have been calculated showing that the re¬ 
versal paths do exhibit a trend towards linear reversal at 
T* but in smaller grains the criticality of the transition 
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FIG. 3. The reversal paths of (a) 2.316 nm, (b) 3.088 nm, (c) 3.860 nm and (d) 4.632 nm systems using a reversing field of 
IT. The reversal paths are calculated at different temperatures slightly below the Curie temperature. For the larger system 
sizes the lower temperatures have not reversed with in the simulation time limit but as the system size decreases at the same 
temperature the system now reverses within the given time, (e) The ellipticity extracted from the reversal paths in (a)-(d) 
using equation 8. The dashed lines show the ellipticity that is expected by the LLB model using the finite size susceptibilities 
shown in figure 2. 
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FIG. 4. The reversal time in a constant 1 Tesla reversing 
field. The solid lines show the analytic minimum pulse dura¬ 
tion which is derived in ref. 2 the parameters for which are 
calculated from our atomistic simulations for the finite sys¬ 
tems while the black line is for bulk. 


is reduced. We find that the expression for the reversal 
time derived from the LLB equation gives switching times 
close to the atomistic results using values appropriate 
for each individual nanoparticle diameter. This suggests 
that a multiscale approach for HAMR media modeling is 
feasible using the LLB equation as a macrospin approxi¬ 
mation to the behavior of nanoparticles, with parameters 
determined for each grain size. Given the importance of 
linear reversal to the HAMR process, this represents an 
important result for simulations of future generations of 
HAMR media as the grain size is inevitably reduced as 


linear densities increase. 
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